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Abstract 

Superkamiokande atmospheric neutrino data for 535 days are analyzed in the 
framework of three flavor oscillations with mass hierarchy. It is shown that 
the best fit point is very close to the pure maximal «-> v T case and Am? ~ 
7 x 10 -3 eV 2 . The allowed region at 90 %CL is given and the implications to 
the long baseline experiments are briefly discussed. 
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Recent data from atmospheric neutrino experiments and especially the Superkamiokande 
experiment [|I],[2j provide very strong evidence for neutrino oscillations. In || and Q the 
atmospheric neutrino data for 414 days |I[ have been analyzed in the framework of two flavor 
oscillations z/ M <-> v r or <r-* i/ s , and it has been shown that both scenarios give a good fit to 
data. In this paper we extend the analysis of [Q] to the case of three flavor oscillations using 
the data for 535 days 0. (For quantitative three flavor analysis of atmospheric neutrinos, 
see [^HlOfl.) In case of general three flavor mixings the data from the CHOOZ experiment 



il| gives a strong constraint to any channel which involves u e , so we include the combined 



X 2 of the reactor experiments CHOOZ, Bugey and Krasnoyarsk JT3 1 in our \ 2 analysis. 

The pattern of mass squared differences with hierarchy can be classified into two cases 
which are depicted in Fig. 1. As in we ignore the smaller mass squared difference, since it 



is expected to be of order 10 eV or 10 eV to account for the solar neutrino deficit [14 



and therefore too small to be relevant to atmospheric neutrinos. Furthermore we put the 
CP violating phase 5 = for simplicity. Then we are left with three parameters (m 2 , #13, #23) 
in the standard parametrization [|T5[] . As has been noted in |7|, two schemes (a) and (b) in 
Fig. 1 are related by exchanging the sign of m 2 , and the case with m 2 > stands for the 
scheme (a) while the one with m? < for (b) in Fig. 1. To evaluate the number of events, 
we have integrated numerically the Schrodinger equation 
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effect in the Earth [0. The way to obtain the numbers of events is exactly the same as in 



|| , and we refer to for details. In || two quantities have been introduced to perform a 
X 2 analysis. One is the double ratio |I7[ 

(iViv e )| osc 

(iV M /iV e )| no _ osc 1 J 

where the quantities N e ^ are the numbers of e-like and /z-like events. The numerator denotes 
numbers with oscillation probability obtained by (|l|), while the denominator the numbers 
expected with oscillations switched off. The other one is the quantity on up-down flux 
asymmetries for a- like (a=e,/x) events (See also |l8|,[19i^Ti Pi .) and is defined by 



a ~ {N-^/N+^)\„ 1 ) 

where N~ v denotes the number of a-like events produced in the detector with zenith angle 
cos < — rj, while N+ v denotes the analogous quantity for cos 6 > r], where rj is defined to 
be positive. Superkamiokande divides the (—1, +1) interval in cosG into five equal bins, so 
we choose rj = 0.2 in order to use all the data in the other four bins. Thus \ 2 for atmospheric 
neutrinos is defined by 



Xatm ^ t 
E 



' R SK_ Rth y , Y SK_ Y th ^ , Y SK_ Y th \* 

6R SK J + { 5Y* K ) + { 5Y e SK ) 



(5) 



where the sum is over the sub-GeV and multi-GeV cases, the measured Superkamiokande 
values and errors are denoted by the superscript "SK" and the theoretical predictions for the 
quantities are labeled by "th" . In |J a x 2 analysis has been performed using the quantities 
R and Y's, or using Y's only. Throughout this paper we use the quantities R and Y's to get 
narrower allowed regions for the parameters. We have to incorporate also the results of the 
reactor experiments. We define the following x 2 : 

CHOOZ / _ \2 Bugey / \2 Krasnoyarsk / \ 2 

j=l,12 \ 3 / i=l,60 \ Udy 3 / j=l,8 \ 3 I 

where X{ are experimental values and yi are the corresponding theoretical predictions, and 
the sum is over 12, 60, 8 energy bins of data of CHOOZ |TT[, Bugey |TjJ and Krasnoyarsk ||13|| , 
respectively. There are 6 atmospheric and 80 reactor pieces of data in x 2 = xitm + Xreactor 
and 3 adjustable parameters, m 2 , 6*13 and 6*23, leaving 83 degrees of freedom. 



Using the same parametrization as that in J7|, the results for the allowed region of the 
mixing angles (#13, (9 23 ) are given for various values of m 2 in Figs. 2 and 3. The results for 
m 2 > and m 2 < are almost the same. It is remarkable that, unlike in the case IF1 of the 



Kamiokande data |T7|], the Superkamiokande data strongly favor <-> v T oscillations. This 
is not only because we have included the combined Xreactor °f the reactor experiments but 
also because the Superkamiokande data themselves favor <-> v T [22| [23 . 



The best fit is obtained for (m 2 , tan 2 13 , tan 2 #23, X 2 ) = (7 x 10" 3 eV 2 , 1.0 x 10~ 2 , 1.6, 
72.8) for m 2 > and (-7 x 1(T 3 eV 2 , 1.0 x 10~ 2 , 1.6, 72.7) for m 2 < 0, respectively. 
Xmin = (Xatm)min + (xLctor)mm = 5 - 3 + 67 - 4 = 72 - 7 indicates that a fit to data is good for 83 
degrees of freedom at the best fit point. The allowed region for m? with # 13 , #23 unconstrained 

is given in Fig. 4, where A X 2 = xl tm + xLctor ~ (xLm + X r 2 e actor)mm < 3 - 5 > 6 - 3 > 1L5 
corresponds to la, 90 % CL and 99 % CL, respectively. The allowed region for \m 2 \ at 99% 
CL is 3 x 10~ 4 eV 2 < |m 2 | < 1.8 x 10 -2 eV 2 . It should be noted that the large m 2 limit is 
excluded because we have postulated the constraint of the reactor data [f20[ . 

Finally, let us discuss briefly the implications of the present analysis to the long baseline 



experiments pi|-|2"6]]. One of the interesting questions in these long baseline experiments is 
whether v e can be observed from <-> v e oscillations which could be present as a fraction 
of the full three flavor oscillations. The probability Piy^ <-> u e ) in our scheme is given by 

P{y, - u e ) = 4|f/ e3 | 2 |^ 3 rsin 2 , (7) 

where L stands for the path length of neutrinos. The factor 4|L/ e3 | 2 |?7 At3 | 2 corresponds to 
sin 29 in the two flavor framework, so by substituting this quantity in a (Am 2 , sin 2 20) 
plot we can examine the possibility of observing v e . The maximum values of the factor 
4|f^e3| 2 |^3| 2 in the allowed region at 90 % CL and 99 % CL are given in Table. 1, respectively. 



In general it is difficult for the long baseline experiments [p4|-[26| to see appearance of u e . 
In particular, for the K2K experiment, which could probe the region of \m 2 \ as low as 



3 x 10 3 eV 2 for «-> v e oscillations [2~7|, it seems very difficult to observe appearance of 
v e and disappearance of <-> has to be searched for at least in the first stage of their 



experiment. On the other hand, if 5 x 10 4 eV 2 < \m 2 \ < 1.0 x 10 3 eV 2 , there is a chance 



for KamLAND [[28|] to see a positive signal in a disappearance experiment of v e <-> u e . 

In conclusion, we have analyzed the Superkamiokande atmospheric neutrino data in the 
framework of the three flavor oscillations with mass hierarchical ansatz. We have given a 
allowed region at a certain confidence level for the mass squared difference and the mixing 
angles. The data strongly favor <-> v T oscillations and therefore the most promising way 
in the long baseline experiments is to search for appearance of — *> v T if v T can be produced, 
or to look for disappearance of «-> if v T cannot be produced. 

The author would like to thank H. Minakata for discussions, R. Foot for many useful 
communications and Center for Theoretical Physics, Yale University for their hospitality 
during part of this work. This research was supported in part by a Grant-in-Aid for Scien- 
tific Research of the Ministry of Education, Science and Culture, #09045036, #10140221, 
#10640280. 
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Figures 

Fig.l The hierarchical neutrino mass squared differences. The scenarios (a) and (b) are 
related by exchanging m 2 <-> — m 2 . They are equivalent in vacuum but physically 
inequivalent in matter. 

Fig. 2 Three flavor analysis of Superkamiokande atmospheric neutrino data and the reactor 
experiments, CHOOZ, Bugey and Krasnoyarsk. Scenario (a) in Fig. 1 is assumed. 
The solid, dashed, dotted lines represent 68 % CL, 90 % CL, 99 % CL, respectively 
for degree of freedom = 3. The right side of each panel corresponds asymptotically to 
pure ffj, <-> u e oscillations and the lower side to pure <-> u T oscillations 

Fig. 3 As in Fig. 2, but the scenario (b) in Fig. 1 is assumed. 

Fig.4 Value Of A X 2 = Xatm + Xr'eactor ~ (Xatm + XLctorLin = Xatm + xLctor ~ 72.7. The Solid, 

dash-dotted, dotted lines represent the scenarios (a), (b) and the two flavor case with 
maximal <-> v T mixing, respectively. 
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Table. 1 Maximum values of the coefficient in the probability 
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Fig.2 
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Fig.3 
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